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The most general class of Bell’s inequalities for A'^-partite d-chotomic systems, with k number 
of obervables, is derived using the formalism of generic correlation functions, which are obtained 
as a fourier transformation of probability spectrums. It is found that they provide full description 
of the probability space together with corresponding operator decomposition for experimentally 
accessible local measurement settings. It turns out that the conditions for the local hidden variable 
(LHV) model are violated by quantum correlated states in general and they constitute the general 
correlation space for Bell theorem. It provides the most general constraints for the LHV model and 
the conditions for the quantum violation. Quantum violation is evident by maximally entangled 
states. 
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Introduction- Since Bell formulated an experimentally 
testable condition for a local realistic model of a di¬ 
chotomic bipartite system [Ij , generalization of the theo¬ 
rem to an arbitrary size quantum system has become one 
of the most challenging topics in the study of quantum 
information science (2 h 13I| . The problem is closely related 
to the possible characterization of macroscopic quantum 
systems quantum key distribution schemes Q, net¬ 
work characterizatio n 1141 a nd entanglement detection in 
many-body systems Mm- 


Conventionally, the generalization of a local hidden 
variable (LHV) model was addressed in terms of “all the 
Bell inequalities” by Peres about two decades ago 0. 
The problem has been challenged in various ways but 
only partial solutions are available so far d, d| ■ Partic- 
ulary, the problem is important because identifying the 
most general class of Bell’s inequalities is related with 
characterization of entanglement and quantum correla¬ 
tions [lB |. However, it is usually known that the prob¬ 
lem is highly non-trivial and is classified as an NP-hard 
problem since it entails the exponentially large parameter 
space of the probability of events. 


In order to detect genuine multipartite entanglement, 
a family of Bell-type inequalities for iV-qudit systems 
were derived ll l| . The main idea is from Svetlichny’s N- 
separability Q and obtained a set of Bell inequalities up 
to a limited dimension, exploiting exhaustive numerical 
approaches [l^- In spite of this progress, the full gener¬ 
alization of the Bell function has not been made yet and 
it is unclear whether the generalization is possible using 
the previous approaches. 


In this Letter, we provide a general method for the 
analytic construction of Bell correlation functions and 
prove that the inequalities are violated by non-trivial 
quantum entangled states. In this construction, it is 
demonstrated that well-known Bell-type inequalities can 
be derived from the generic form of the correlations. Con¬ 
sequently, the analytical condition of the LHV model is 
provided. 


General class of Bell’s inequality - The class of Bell 
inequalities is determined by the three parameters N, 
k and d; N local parties measure their systems with k 
possible choices of observable which result in d differ¬ 
ent measurement outcomes respectively. Once N, k and 
d are dertermined, we can write conditional probabili¬ 
ties which constitute information about the system [T^: 
a set of functions, which constitute a catalogue of full 
information, take the form of conditional probabilities 
p{ai,a 2 , • • • , aAr|mi, TO 2 , • • • , itin) with I < < d and 

0 < ruj < fc — 1, where aj is an integer number indicating 
the measurement outcome index for a particular choice 
of the observable mj. A Bell correlation function, which 
may be measurable experimentally, is derived from the 
one-to-one correspondence between the probabilities and 
multi-partite high-order correlation functions. With the 
Bell functions, we can then formulate generalized Bell 
inequalities. 

Mapping the measurement outcomes to d different val¬ 
ues, the most general correlation function can be written 
as 


Eni {mi,m 2 , ■■■) 



( 1 ) 


where 1 < j < N is the site index. We define a shorthand 
notation En,^^n 2 ,---{’n^i,tn 2 ,- ■ ■) '■= Efi{rn). By mapping 
the measurement outcomes Aj{mj) to one of the complex 
values among the d root of unity, Aj{mj) = uj°‘Araj)^ 
correlation function takes a complex valued number with 
w = exp(i27r/d) and aj € {1,2, ••• ,d}. The function 
allowing predetermined values of a measurement, is eval¬ 
uated by the expectation of the measured values (out¬ 
comes) 


d 

En{m) = ^ a;”'“p(d|TO) (2) 

{“31 = 1 

where a. := (q:i(toi), 02 (^ 2 ), •••, Q;Ar(mfe)). aj{mj) 
is determined by the physical process of measurement 
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whose full structure is hidden and it is usually encapsu¬ 
lated in the form of abstracted variables, conventionally 
denoted by A. 

It is notable that the correlation that is to be obtained 
from the probabilities can be defined differently in gen¬ 
eral. That is because the combination of measured values 
can be arranged in an arbitrary manner as there is no 
preferred order of indices for the outcome sequences. In 
order to consider the full sequences of arbitrary combina¬ 
tions, it is necessary to introduce extra integer index c for 
a type of correlation: Efi^{rh) = X]{aj}=i 
where ftc = (cini, C 2 n. 2 , 0371 , 3 , •••) with Cj G { 1 ,- 1 }. 
The subscript c takes an integer value between I and 
2^ and specifies the type of possible correlations. Sub¬ 
sequently, it is possible to recover all the spectra of rel¬ 
evant probabilities p(a| 77 i) = 


where the total number of distinguishable probabilities 
are d^k^. The equality provides the one-to-one corre¬ 
spondence between probabilities for measurement out¬ 
comes and the high-order correlations. Without any 
other constraints, the probability function satisfies (i) 
positivity p(a| 777 ,) > 0 , Vaj,V? 7 ij and (ii) the normaliza¬ 
tion condition X]{o! } p(<5|w) = I as the correlation func¬ 
tion is required to satisfy \Efi^{m)\ < I and ifj- {m) = I 
for any choices of m with dc = (cid, C 2 d, ■ ■ ■)■ 

If we generalize the interdependency of the probabili¬ 
ties from different measurements, due to incompatibil¬ 
ities, it is possible to derive a Bell-type inequality in 
the most general situation. This is main result of this 
work. Using the definition, the most general Bell func¬ 
tion is given by 


G C 

N,k,d 


d k—1 


d-1 




where c.c. denotes complex conjugate. The coefficient 
gg is a real function governing the linear sum of the 
probabilities and f{nc) is a complex weighting function 
for the correlation. For the LHV model, the average is 
taken over the hidden variable A for the measurement 
which determines the value of the measurement function 

Aj. 

The choice of the function /{fie) = f{cini,C 2 n 2 , ■ ■ ■) 
specifies the weight of each high-order correlation func¬ 
tion in the sum and is related to the coefficient g'b ^ as 


= 2 Re 




(4) 


From its fourier analysis, it is also possible to show that 
the choice of function /(fic) generates all the possible 
combinations of the correlations Efi^ (m). The decompo¬ 
sition of the generation function is obtained when the full 
correlations have been taken into account. Mathemati¬ 
cally, it means that ~ ^ 

when Uj = d, Vj. 

The decomposition of the correlation into the proba¬ 
bilities in Eq. ([3|) provides the most general correlation as 
it constitutes “all the Bell inequalities” Q . The only dif¬ 
ference of this construction from the original formulation 
is the condition of homogeneity. It means that the num¬ 
ber of measurements and the dimension of each party 
are chosen to be identical. It can be shown that all the 
known Bell functions within the homogeneous condition 
can be derived as a special case of the function G in ([3]) . 

For the LHV constraint, the first decomposition in 
using probabilities is directly linked to the modi- 



k-l 

rrij =0 


+ c.c. 


(3) 


fied version of Farkas lemma M- The statement is 
that if J2s,rh9s,rhIlj^jio^j\'mj,X) < Blr is satisfied 
by all A and the local realistic (LR) bound Brr then 


G%,k,d = Ea.m 55 .^P(a|rn) < Blr where p{d\rn) = 
Ea pW rij {aj\mj,X) and p(A) is a positive func¬ 
tion, satisfying Ea P('^) ~ formulation, it is 

straightforward that the LR bound can be obtained by 


the Farkas vector g^ ^ as max^ 


Em 9a ,rn 


There¬ 


fore, the bound for the correlation G5^r k d obtained as 


G C 

N,k,d 


< 


Blr = max 

a 



(5) 


which provides the most general criteria for the prob¬ 
abilities allowed by the LHV model. Evaluation of the 
bound requires the functional optimization over the mea¬ 
sured values, and can be achieved by the specification of 
local parameters under the functional constraints. The 
analytical values can be obtained efficiently if one follows 
the optimal counting method described explicitly in . 

The decomposition of the correlation function for 
possible local measurements under the LHV model is 
nonetheless trivial. It allows quantum characterization 
of the correlation when the coefficients g^ ^ as well as 
/{fie) are appropriately chosen. In the Bell correlation 
function ([3]), the general decomposition of measurements 
with the weighting factor f{nc) is obtained. In the follow¬ 
ing, we demonstrate the derivation of known Bell func¬ 
tions through the specification of g^m and f{nc). The 
condition for quantum violation will follow. 

Derivation of all known Bell inequalities - We first 
show that the function ([3]) reduces to the Clauser-Horne- 
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Shimony-Holt (CHSH)-Bell inequality 0. By taking the 
generation function = (_x)ai+a 2 +mim 2 ^ which 

results in < 2, according to the bound given 

by Eq. ([5]). Here, one can obtain = (1 — 

11 ^ 712 ,i/2- The correlation function becomes Collins- 
Gisin-Linden-Massar-Popescu (CGLMP) function @ for 
(2,2,d)-class system when 


CGLMP 

ya^rh 


d-1 

'y 

k=0 



X S'^(a2 — ai 


2k \ 
d-l) 

— k — z{mim2)) 


( 6 ) 


where S‘^ is kroneker delta function up to modulo d 
and 0 ( 7711 , 7712 ) is a binary function mapping z(0, 0) = 


< max, 

function is given by 

CGMLP _ _ 1/2 


f 

J 1 


0 and 

0(0,1) 

z_-/m ydi,rh 

d-l 



-I] sec 

UTT 

UJ 

I2d\ 


= 2. The weighting 


n—1 


-'ni—nn2——n 


(7) 


whose detailed derivation is shown in M- The quantum 
violation of the LR bound can be found in the original 
works [2, Q and can be recovered in the framework of 
generalized correlation formalism in an analytical way. 

We demonstrate that we can derive Bell func¬ 
tions for multipartite systems such as Mermin and 
Zukowski-Brukner (ZB) functions S 0. The Mer¬ 
min function G/f 2 2 obtained by assigning = 

(—1)^7 “^’Re[(7)^jwhen we have/'^(l) = 1/2. From 
the form of gs,irn the Mermin bound can be found as 
Gn 2 2 — for odd N and 2 2 — for even 

N. Ultimately, the most general (N,2,2)-class ZB corre¬ 
lation can be obtained when f^^in) = (1 — i)/2 as 


G%^-^2= X! P{'mi,m2,-■ ■ ■ ■ E{mi,m2,-■ ■) 

{rrij }—0 

( 8 ) 

where the parity function P(mi,m 2 , ■ ■ ■) € {—1,1} takes 
its value 1 for [{J2j^j) mod 4] S {0,1} and —1 for 
[(Sj-mj) mod 4] G {2,3}. The LR bound of the func¬ 
tion can be found as G^ 2 2 — 2^ from the probability co¬ 
efficient = (—1)^7'^^“^Re[(l — 7 )( 7 )^ 2 through 
the maximization of 9m,s over a proper choice of 
(ai,a 2 ,'--) following its derivation in Eq.(|3]). A more 
general constraint of Bell’s theorem for the class of full 
correlation can be obtained as one combines G^'^^ for 

all the possible c as X]c=i 

For the case of (2, k, 2)-class Bell test scenario, the 
general correlation function has been inspected by Ep- 
ping et. al. [l^ and can take the form — 

Smi,m 2 '^™i.™2-®(mi,77i2) where E{mi,m 2 ) is the first- 
order correlation function. In this case, Pmi,m 2 = 
^ ^ such that the coefficient j3 can be 


related to the probability coefficient as = 

(—l)(“7+“2^/3mi,m2- Thus, the LR bound of the function 
is obtained as 


G 


EKB 
2,k,2 — 


< max 




(9) 


= 2|/(l)|max 


cos(7ra) cos (j)f 


TTm\ 

~ ) 


= 2|/(1)| COS0//sin^(7r/2/c), for 0 < 0/< 7r/2/c 


where d = ai -I- 02 , fh = mi + m 2 and Of is the phase 
factor of the complex function /(I) as /(I) = |/(l)|e®®^. 
The second equation is obtained after optimization by 
counting the largest cosine terms. Quantum maximum 
of G 2 ^ 2 ’^ i® given by the singular value decomposition 
of the ’l3 matrix, fc||/3||2 = fc^|/(l)|, as argued in [i0 . 

Quantum violation - For quantum systems, the mea¬ 
sured variables are expressed by eigenvalues of an op¬ 
erator whose expectation corresponds to the statistical 
average of measurement values. In that case, the decom¬ 
position in Eq.® can be represented by measurement 
operators whose explicit form is given by 


= X]w“|A„(777j))(A„(mj)| (10) 

Q! 


where d-dimensional orthogonal bases \Aa{mj)) are cho¬ 
sen to satisfy {Aa{mj)\Ap{mj)) = 5a,i3- The bases 
\Act{mj)) can be obtained as a linear combination 
of the orthogonal computational bases, \Aa{mj)) = 
1^) where the mj-th basis has been 
obtained by the phase shift of the fourier transformed 
state whose phase shift is distributed from 0 to fc — 1 
evenly. In the sense that the measurement bases are 
evenly distributed in Hilbert space, the bases can con¬ 
stitute a maximal test. From the bases, one can show 
that the spectral sum of measurement operators become 
ladder lowering operators as 

k— 1 

rrij—O j3 

which corresponds to the nj-th power of a lowering oper¬ 
ator Jj = J 2 i 3 |/3)i(/3 + 1| for a high dimensional state as 
{JjY^\a) = \a — Uj). In addition, a phase shift operator 
Pv acting on the orthogonal computational bases gener¬ 
ates an extra phase Pu\a) = and can be used for 

the local unitary transform on the lowering operator as 
P^jJ^Piy = j"7. The phase shift operation is effec¬ 

tive in order to obtain the different measurement bases 
that produce the correlation values beyond LR bounds. 

Together with the sequence of the local phase shift op¬ 
eration P ^., the generalized Bell function for a quantum 
state can be obtained in terms of high-order correlation 
functions. If there is a state whose expectation goes be¬ 
yond the LR bound, the state cannot be described by the 
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LHV model. With the measurements and local rotations, 
the correlation function of a quantum state becomes 

d-l IN \ 

E (12) 

n=l \i=l / 


as long as Blji < Qm- Explicit criteria in a c oup le of 
special cases and their analysis can be found in 

As an example, the quantum upper bound of CGMLP 
equation for ME state can be analytically formulated. 
For ME state, the quantum expection of general Bell 
function becomes 


where v = {vi,V 2 ,‘ ‘') represents the composite compo¬ 
nents of local phase shifts at each site. 

In the following, we show that the LR bound is violated 
by a simple symmetric quantum state with an appropri¬ 
ate choice of /(n). When the powers of the lowering oper¬ 
ators at each site are uniform, ni = n 2 = ■ ■ ■ = n, the N- 
partite maximally entangled (ME) pure quantum state 
IV') = |q,)®./v gjygg quantum correlation 

(1 - ^) + c.c. with an 

arbitrary local parameter Vj representing a choice of mea¬ 
surements. If p is chosen to satisfy that —n pj = 
Arg[/(n)], the value of the quantum correlation is upper 
bounded by the quantum maximum for the ME state 

GN^,d ^ X;Ei (1 - l) l/(^)l = Qm- Thus, an ap¬ 
propriate specification of f{n) results in the violation of 
LR bound 


d-l 

= |/(n)| (13) 

n—1 

the bound Blu is also determined by the function /(n), 
as shown in Therefore, the general constraint for 
Blr < Qm can be obtained from the appropriate choice 
of the weighting factors. If |/(n)| is either a constant 
function or a monotonically increasing function with re¬ 
spect to n/d, then a general trend of violation, Blr < 
Qm, can be obtained. Generally speaking, the values of 
Blr determined from f{n) provide the criteria for the 
LHV model and constitutes the generalized Bell function 


o , nCGMLP 


d-l 


d — 1 


E (i - 0) 


n—1 


'nn' 

.M. 


(14) 


for any value of d. Thus, the state is non-local although 
the value is not quantum maximum. It is known that the 
maximal violation of the CGMLP inequality is obtained 
by partially entangled state. Indeed the maximum viola¬ 
tion value can also be derived from the current formalism. 

Remarks- In this work, we derived general criteria 
of the LHV model analytically and provided examples 
that violate the condition. The derivation was done 
through the one-to-one correspondence between the gen¬ 
eral probability space and the correlations of arbitrary 
high-order moments. We proved that the decomposition 
of the correlation function using joint probabilities gives 
the straightforward LR bound. The decomposition by 
full correlations of high-order moments can be used for 
non-trivial quantum violations under the provided set¬ 
tings of measurements. The result sheds the light on the 
general characterization of quantum correlation in an ar¬ 
bitrary number of high dimensional systems. 
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